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TWO-COMPONENT INTEGRABLE GENERALIZATIONS OF BURGERS 
EQUATIONS WITH NONDIAGONAL LINEARITY 

DARYOUSH TALATI AND REFIK TURHAN 


Abstract. Two-component second and third-order Burgers type systems with nondiagonal con¬ 
stant matrix of leading order terms are classified for higher symmetries. New symmetry integrable 
systems with their master symmetries are obtained. Some third order systems are observed to 
possess conservation laws. Bi-Poisson structures of systems possessing conservation laws are given. 


1. Introduction 


Systematic classification of multi-component integrable equations is initiated by Mikhailov, Sha- 
bat and Yamilov in [T]. They completely classified second order, two-component systems of form 

^ = A{u)u^^ + F{u,u^), det(A(u))/0, 

possessing infinitely many conservation laws. One of their extensive results was the fact that all 
the integrable cases of the class considered were those systems that can be written by generalized 
contact transformations in the form in which the coefficient matrix of leading order terms A = J 2 
where is the second of the three constant matrices 


Ji = 


1 

0 



1 0 
0 A 


A/1; J 


(a,e) 

3 


a e 
0 a 


,e^0. 


Up to an overall factor which can be absorbed by the evolution parameter t, and permutation of 
rows which corresponds to rewriting the system with u and v interchanged, the above matrices, 
which are deduced from the canonical 2x2 Jordan matrices, exhaust all possibilities for the 
constant coefficient matrices of linear (in leading x-derivative order M) terms of any quasilinear 
two-component system 

— = Aum + lower order terms, 
dt 

under the linear change of dependent variables (similarity transformations of constant matrix A). 
Also, up to a rescaling of t, a = 0 and a = 1 suffice for systems with A = and nonzero e 

can be adjusted to any desired nonzero value by a rescaling of v. So for quasilinear systems with 
nondiagonal linearity, the constant matrix A can be chosen, without loss of generality, to be 
and to represent the nodegenerate and the degenerate cases respectively. 

Multi-component generalizations of second order 


du 

dt 


B2 [ r ] — Uxx ^UUx ? 


and third order 


R.3[r] — Uxxx 




scalar Burgers equations have been the subject of various systematic integrability classifications. 
Among these, those rewieved in [2] and the two-component cases of Svinolupov [3] were classifi¬ 
cations with assumption of matrix A = Ji. In the classifications by Sanders and Wang [1], and 
two-component cases of Tsuchida and Wolf [5] the case A = with arbitrary nonzero A are 
classified. Burgers type systems classified in [S] were systems with matrix A similar to Ji and 
Here we mention papers pertaining to multi-component generalizations of Burgers equation only. 
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For the other integrable systems and their properties, we refer to review papers mmu and the 
references therein. 

So far the only known integrable system with nondiagonal A is the Karasu(Kalkanh) system 

— Wxxx Zxxx “b 2WxZ + 2wZx , 

\ — Zxxx OwWx “b GzWx “b 3wZx “b 2zZx , 

obtained in Painlave classification [9]. A recursion operator is given for the system in [lOj and its 
bi-Poisson formulation is given by Sergyeyev in m- System ([T]) is related to the system 


f — Uxxx “b 3uUxx “b 3u Ux “b 3u^ Vxx “b 2vUx , 

\ — '^xxx 3uVxx “b QvUxx “b 3'iix'^x “b Vx 2vVx , 


by the noninvertible Miura transformation w = Ux, z = v + The reduction u = 0 of the 

transformed system m is nothing but the scalar third-order Burgers equation ^ = B^lu]. 

As a general rule, as the number of components and/or the order of a system increase, the 
computations needed for a higher symmetry rapidly goes beyond the available computing power. 
In such a situation it is not possible to consider a class of systems in its full generality. One of 
fruitful approaches to specify a rather restricted class of systems is to consider systems made up of 
only polynomial terms selected according to a certain weighting scheme. This approach is based on 
the observation that the vast majority of known integrable equations can be written in polynomial 
form in which they are homogeneous in a certain weighting scheme defined by invariance under 
scaling 

(x, t, u, v) {a~^x, a~^t, a^^u, a^^v), a > 0 (g M). 

Such a system is said to be (Ai, A 2 )-homogeneous of weight fi. In this terminology, system ([T]) is 
(2,2)-homogeneous of weight 3 and its transformed form ([2]) is (1,2)-homogeneous of weight 3. 

Here we classify (1, l)-homogeneous systems which are left unclassified in terms of the matrix of 
leading order terms in HEIE]. I.e. we classify (1, l)-homogeneous of weight 2 and 3 class of (Bugers 
type) systems having nondiagonal constant coefficient matrix of leading order terms A = to 

admit a symmetry from the class of (1, l)-homogeneous of weight 1 and/or 2 higher than that of 
the class itself. We give complete list of systems having nondegenerate i.e. the case with 

a 7 ^ 0 and discuss some degenerate (a = 0) cases that arise as certain limits of nondegenerate cases 
obtained. 

Two systems of A-component evolutionary partial differential equations, or shortly systems, 
with characteristis P and Q 


(3) S = PM . 

are said to be compatible if their flows commute 


( 4 ) 


d du _ 

dt dr 


du 

dr 


Q[u] 


d du 
dr dt 


Here, u = (n^, • • ■ , and the dependent variables n*, i = 1, 2,3, • • ■ ,N are assumed to be 

dependent on independent variable x and evolution parameter t. P[u] means that the components 
P* of P are differential functions depending smoothly on the independent “spatial” variable x, 
dependent variables and finitely many x-derivatives of tt® which will be denoted as u\ = 

For small values of /c < 3, we often use instead of u\,U 2 ,u\. The highest order of 

x-derivative of tt* in a differential function is referred as the order of the differential function. 

Using chain rule in the left hand side of compatibility condition (j3|) one gets 


d 

dt 


Q = 


E 



ApQ 


where 


E <+i 

l<i<N, 

A:GZ>o 


d 

dul ’ 
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is the total x-derivative and 


Xp 




d 



is the evolutionary vector field with characteristic P. Evolutionary vector fields with at most first 
order characteristics generate point transformations. Those with higher order characteristics are 
referred as generalized evolutionary vector fields. 

By a similar identification in the right hand side, compatibility condition Q reads 


= Is- = »• 

Therefore evolutionary systems in (l3|) are compatible if (and only if [12]) the Lie bracket ^[Xp,Aq] 
of their corresponding evolutionary vector fields vanishes. An evolutionary vector field commuting 
with the characteristic of a system is defined as a symmetry of the system. If a (dispersive) system 
possess infinitely many symmetries (of arbitrarily higher orders) then the system is defined to be 
symmetry integrable. A collection of mutually commuting evolutionary vector fields constitute an 
abelian Lie algebra which is also referred as symmetry algebra. 

A particular system which is determined to possess a higher order symmetry is a natural subject 
of further investigations for structures like recursion operator |12] or master symmetry mm 
existence of which indicate existence of infinitely many higher symmetries. Because by definition a 
recursion operator maps a symmetry to another symmetry of a system and master symmetries do 
the same by their adjoint action on a given symmetry. 

Some symmetry integrable systems can be included in an infinite hierarchy of compatible systems 
each of which is a conservation law 

^ = Pn[u] = nSjhn+i[u]dx = }CSjhn[u\dx, n = -1,0,1,2, • • • , 

CLtn 

written by two compatible Poisson structues % and /C, each, and sum of which are skew-adjoint oper¬ 
ators satisfying Jacobi identity. Here, <5 denotes the variational derivative with components {6)i = 

Hamiltonian functionals in involution = 

f 6hn ■ 'Hdhmdx = 0 = f 6hn ■ tCdhmdx = {Jhn, f hm}K with respect to both Poisson brackets corre¬ 
sponding to each Poisson structure. Construction of such infinite hierarchy of compatible systems 
and the corresponding Hamiltonian functionals in involution, i.e. the Lenard-Magri scheme, in 
the case of local Poisson structures is rigorously addressed in [15] for single-sided and in [MKH] 
for multi-sided Lenard-Magri schemes. The case of nonlocal Poisson structures is treated in |18j . 
where a completely algorithmic criteria (i.e. without needing any user intervention) is given for 
checking whether a nonlocal skew-adjoint operator /C which is a ratio of local operators A and B 
as 1C = A{B)~^ satisfies the Jacobi identity or not. We do not go any further details of this broad 
subject here but refer to the mentioned references and to the book |19| by Dorfman for a survey of 
early results. 

We present resnlts of onr computer classification of (1, l)-homogeneous of weight 2 and 3 sys¬ 
tems with nondiagonal linear part admitting a certain type of higher symmetry. We give master 
symmetries of systems determined to have only higher symmetries, and bi-Poisson formulation of 
those that turned out to possess conservation laws also. 


2. Second-order systems 

The class of two-component (1,1)-homogeneous of weight 2 systems with undetermined constant 
coefficients c* have the form 

^ = c\uxx + clvxx + cluUx + cluVx + CgUtta, CqVVx + 

-|- -|- Cgun^ + , 

^ = cjuxx + cjvxx + cluUx + c\uVx + clvUx + C^Wx + 

+ clu'^V + C^uv"^ + . 


( 5 ) 
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The class of (l,l)-homogeneous of weight 3 systems with undetermined constant coefficients Ij is 
^\Uxxx T J‘2^xxx ^^UUxx T li^UVxx T T ^qI^Vxx T ^J^x ^g^x^x 

+ Igvl + ligu'^Ux + liiU^Vx + I12VUUX + ll^^UVVx + ij^v'^Ux 
+ ll^v'^Vx + + lljU^V + llgU^v‘^ + l\gUV^ + , 

^l^xxx T ^2^xxx T IgUUxx T l^UVxx T ^^VUxx T ^Q^Vxx T ^j'^x ^g^x^x 

+ IgVx + tlgU^Ux + ihu^Vx + I 12 VUUX + if^UVVx + tl^V^Ux 

+ if^v'^Vx + IiqU'^ + if'fU^V + + ifgUV^ + ^20^^ • 

The class of (l,l)-homogeneous of weight 5 systems contains about 60 terms in each component. 
Therefore we omit writing this class explicitly. 

Imposing compatibility condition (SD among the classes of systems m and dS]), we obtain a 
system of (initially bilinear) algebraic equations among the undetermined constants c*- and /*’s. 
Each solution of the system of constraints on constants c)- and Ij determines a system and its 
higher symmetry. The system of constraints are solved by the CRACK package [20]. Among the 
systems obtained to possess a higher symmetry, the uncoupled ones and the triangular ones, i.e. 
those that reduce to successive scalar equations, are discarded as they are trivial. The remaining 
systems are searched for conservation laws by the package m- 

Solutions of the compatibility condition (Hj) among the classes (0) and dHj) with cj = c| = 1, 
cf = 0, C 2 (= e) / 0, (i.e. for A = is given in the following proposition. 

Proposition 1. Any (1,1)-homogeneous of weight 2 system of form ([5|) with nondegenerate con¬ 
stant coefficient matrix of leading order terms having one dimensional eigenspace, possessing a 
symmetry from the (1,1)-homogeneous of weight 3 class of systems dGj) is equivalent, by a linear 
change of variables u and v and rescaling of x and t, to the following equation 

(7) / ft=B2[u] + e{vx + vu)^, 

I ^ =Vxx + 2,uVx - e{uv‘^ + vvx ), 

with e 7 ^ 0 . 


( 6 ) 


^ du 
dr 


dv 

dr 


Remark 1. The characteristic of system ([7|) appears like a linear combination of the characteristics 
of the triangular system 


( 8 ) 


/ W = B2[u], 

I ^ = Vxx + 2uu 




and the characteristic of the system 


(9) 


/ M — T VUx T UVx 
I ^ = -UV^-VVx, 


which has the degenerate nondiagonal matrix A = Indeed, systems © and are 

all compatible with each other. However, despite the fact that these three systems have common 
symmetries (of higher orders), each of these systems individually have different symmetry algebras: 
Some elements of the symmetry algebra of system © are not admitted by systems dHj) and © as 
a symmetry. Similarly, system dSj) (and ([9 |) ) has its own symmetries not admitted by the the other 
two. Therefore, the eonstant e appearing like an arbitrary constant in a linear combination in the 
characteristic of system © cannot be regarded as sueh, and e 7 ^ 0 is not a removable artifaet of 
the formulation of the classification problem. 


Remark 2. The Hopf-Cole transformation u = ^ and v = ^ not only linearizes system © but 
makes it a triangular system. 


Remark 3. System ([ 8 |) becomes a first order equation of hydrodynamic type by u = Ux. 
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3. Third-order systems 

Setting l\ = I 2 = 1) = 0, l 2 {= e) 7 ^ 0, (i.e. A = in ([ 6 ]), and imposing compatibility 

condition dH among the class of systems ([ 6 |) and the class of (l,l)-homogeneous of weight 5 systems, 
we obtain the systems listed in the following proposition, as solutions of conditions imposed on the 
undetermined constants. 

Proposition 2. Any (1,1)-homogeneous of weight 3 system of form ([ 6 |) with a nondegenerate 
constant coefficient matrix of leading order terms having one dimensional eigenspace, possessing a 
symmetry from the (1,1)-homogeneous of weight 5 class of systems is equivalent, by a linear change 
of variables u and v and rescaling of x and t, to one of the following eight equations: 


( 10 ) 

( 11 ) 


/ W = ^3[u] + e{vxx + ‘2uVx + VUx + vfv)^ , 

1 ^ = Vxxx + 3UxVx - SvUUx - e{vvxx -Vx + V^Ux) , 

I w = ^3[u] + e{vxx + ‘^UVx) ^ , 

\ — Vxxx i^uVxx T Vx T ; 


( 12 ) 

(13) 

(14) 

(15) 

(16) 
(17) 


/ ^ + ‘i^{vxx + “iuvx + 2vUx + 2u^v)^ , 

1 ^ = Vxxx + ^UVxx + GUxVx + 2>U^Vx “ C^'ivVxx - + 8v^Ux + SuVVx + , 

/ ^ = B3[u]-\-e{vxx + 2uVx + vUx + u‘^v)^, 

\ — Vxxx T QVjxVx (^{yVxx Vx; T V Uxf) , 

I ^ + ‘ie[vxx - vux)x, + 2>e^[^vVxx + 8vl - v'^Ux) + 12e^v'^Vx , 

1 ^ = Vxxx - 8{uVx - u^v)x + 3e{8vVxx + 3vl - 2{v‘^u)x) + 21e‘^v^Vx , 

f ~dt ~ T ^{vxx T VVx) x- 1 
\ — Vxxx 3uVxx T 3vxUx T 3u Vx T 2eVx ., 

/ W = ^3[u] + e[vxx - uvx - 2vUx - 2u‘^v)^ , 

1 ^ = Vxxx - 3{uVx - u^v)x + e(2vvxx + vl- GuVVx - 4v^Ux) , 

/ ^ = -^sN + 3e{2vxx - 3vUx - 3e^{y‘^u)^ , 

1 ^ = Vxxx - 3{uVx - u'^v)x + 3e{3vvx - + 3e^v‘^Vx , 


where e 7 ^ 0 . 


None of the systems given in the above proposition is in the symmetry algebra of a second order 
system. Third order systems which are identified to be a symmetry of a second order system are 
discarded. 

All systems (nnn-diTi) reduce to scalar third-order Burgers equation ^ = B^lu], by u = 0 
reduction. 

Linear transformations of dependent variables preserve homogeneity of the classes considered 
but they are too restricted as an equivalence criterion. Therefore, hereafter we investigate each of 
the systems in the above proposition for a possible transformation, not necessarily invertible that 
can relate the system to a trivial one, to a known system, or to a nicer form of it; for a master 
symmetry, and for bi-Poisson formulation. 

3.1. System (fTOjl . System (fTOl) transformed hy w = and . 2 ; = e~ ^{ve ^or equiva¬ 

lently in three stages u = ^,v = ^, Sx = wz, becomes the triangular equation 

/ 717 — vJxxx T ^(vjz)x ) 

\ — Zxxx T (^^Zx , 


which is trivial. 
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3.2. System (fTTT) . System m transformed by u = ^ and Vx = pq, becomes 

f ^ = Pxxx + e{p‘^qx + 3qppx), 

\ — qxxx Px ^pqqx} ; 

which is the Sasa-Satsuma system [22]. Poisson-symplectic formulation of the Sasa-Satsuma system 
is given by Sergyeyev and Demskoi in [23| with some misprints and by Wang in [24] . 

3.3. System (|12p . System ()12p has the master symmetry 

]YjiIl2]i ^ f + 2uxx + 5uux + u^ + e{3vxx + Qvux + 8uvx + iu'^v) 

V + Vxx + 3vux - 2e{vvx + 2v‘^u) 

So one can generate symmetries of arbitrarily high orders by using ad^jrg successively. Therefore 
system (d, obtained by acting on x-translation symmetry, is symmetry integrable. 

Similar to the case in Remark [H the coefficients of and in the characteristic of system m 
and the characteristic of system m as a whole, are all symmetries of each other. But their 
symmetry algebras do not overlap. 



3.4. System (fTHI) . System m is a symmetry integrable system having the master symmetry 
^ + 2uxx + ^uux + + \e{‘^Vxx + ^vux + 7uvx + 3u^v) \ 

V x^ + Vxx + ‘2^uvx - \e{yvx + uv^) ) ' 

Acted on x-translation symmetry, gives system m- 


3.5. System (fl^ . Symmetry integrable system m is obtained by acting ctd-^ to ^-translation 
symmetry where the master symmetry is 


M® 


, du 
" dt 


+ 


X + Quux + 

x'ddL _l_ ^.71 _ 

X dt ^ 2 


+ e(6ua;a; — Avux — uVx) -|- e^{l8vVx — v‘^u) + 
3uvx -I- 3u^v -I- e(l3OTa; — 6u^n) -|- 7e^v^ 


3.6. System (I15p . Transforming system (1151) by ^ for convenience and setting without loss 
of generality e = 1, we obtain 


(18) 


^ = Bslu] + {Zx -t uz)x , 

— Zxxx T 3 (^UxZ UZx 'll z'jx AiZZx^ 


for which we have the following proposition. 


Proposition 3. System ()18ll is bi-Poisson 


^ ~ T~l.dJz{ux + ^u^-\-^z)dx = JCSf^dx, 

with the compatible pair of Poisson structures 

( -\D D‘^ + 2Du\ ( /Cii /C 12 \ 

V + 2zx)’ V -^*12 ^22 j ’ 

where 

/Cii = —D^ + u^D -|- uux — UxD~^Ux, 

/C 12 = + AuD^ + {9ux -|- 3z -|- 5u^)D^ + {7uxx + 7>Zx + Wuux -|- 3uz + 2u^^D 

T (f^'ii'xxx T ^‘ii'i^xx T ^^x A~ 6rx Ux -1- ‘^Zxx T 3zux -1- 2uZxJ 'iixD Zx-, 

/C 22 = QzD^ A- dZxD"^ -|- [7zxx + l2uxZ — 12zxU -|- 12u‘^z -|- 9z^)D 

T (f^Zxxx ^'ll'Zxx T ^ZUxx T ^2uUxZ Zx “i” 9zZxJ ZxD Zx- 

Here, D~^ denotes the formal inverse of D and ICI 2 is the formal adjoint 0 //C 12 . 
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Skew-adjoint local operator T-L satisfies Jacobi identity. Therefore it is a Poisson structure |15j . 
To show that the skew-adjoint nonlocal operator /C is a Poisson structure and compatible with T-L, 
it suffices to verify that the fractional decomposition of nonlocal operator aT-L + K, = AB~^, by 
local operators A = {oiH + JC)B and 

1 


B = 


-D 0 


0 




satisfies the Jacobi identity in terms of operators A and B in [18] regardless of the value of constant 
a. The case a = 0 corresponds to Jacobi identity for /C and a 7 ^ 0 to the compatibility of % and 
1C. These completely algorithmic computations which does not require any user intervention, are 
carried out by computer. 

Recall that the case of arbitrary nonzero e can be recovered simply by transforming all above 
quantities by z —>■ e^. 

3.7. System (fT 6 |) . As in the previous system, for the system (11611 with e = 1, we have the following 
proposition. 

Proposition 4. System ()16p is bi-Poisson 

— = Pbf ( — UxVx + Svuux — ^v'^Ux — ^Vx + vv? — dx = 1C5f uv dx 

with the compatible pair of Poisson structures 

D — |JJ o v{D -|- 2u)~^ 

L> + |(T) - 2u)-i o vD -|((T> - 2u)-^v^ + v^{D + 2u)-^) 

and 1C having components 

Kii = — u^D — uux J- UxD~^Ux, 

IC 12 = + 2{u — v)D‘^ -|- {3ux — 3vx — 2uv -|- u^)D 

+ {Uxx — Vxx + ‘2.UUx — 2vUx — UVx) + UxD~^Vx, 

V' V'* 

^21 — -^125 

JC 22 = —D{vx — 2,uv -|- 2v'^) — {vx — 2uv -|- 2v‘^)D -(- VxD~^Vx. 

A fractional decomposition 1C = A]c{B]c)~^ of nonlocal skew-adjoint operator 1C is obtained by 
the local operators 


H = 


Bk = 


^-D 0 

0 


and Afc = ICBjc- Skew-adjoint operator 1C satisfies the Jacobi identity in terms of the local operators 
Aic and Bjc in [T8|. Therefore /C is a Poisson structure. A fractional decomposition for the skew- 
adjoint operator P is obtained as P = A'u{By)~^ by the local operators 

\D o vR{D - 2u) {D^ + 2Do{u- |)) o R{D - 2u) 

{D + |u) ovR{D- 2u) - I -|(u2R(T) - 2u) - l) 


Au = 


Bh=( 


0 


0 {D + 2m) o R[D — 2m) 

where R = ' Local operators A-h and By^ satisfies Jacobi identity for P [18]. Therefore P 

is a Poisson structure. For the compatibility of P and K. one needs to verify the Jacobi identity for 
P-\-lC. For P-\-lCwe. have a fractional decomposition P-\-lC = Ay+i<z{Byy]c)~^ by local operators 
Ay+K. = {Ay -I- lCBy)S, By+y = ByS with 

S = ( ^ 

V 0 io 

where F = {vUxR)^ + 2vuUxR and G = {vxxR)x ~ ‘2vx{uR)x — ^u'^VxR- The Jacobi identity for 
P + 1C, m terms of Ay+y and By+y is verified to be satisfied. So Poisson structures P and 1C are 
compatible. 
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System (fTHj) . transformd hy u = and v = WxZx becomes 


(H — Wxxx “ 1 “ (-WxiWxZxx WxxZx') 
dz 


dl — Zxxx ^ZxiWxZxx 'UJxxZx') : 
which can be complexified through ip = w + iz as 

= P^xxx + ~ '4’x'lpxx)- 

This equation is the potential form of an equation given with its Lax pairs by Tsuchida in 

3.8. System (fT7|l . For the convenience of Poisson structures, besides taking e = 1 in system (dD, 
we transform it by re = + u), z = — u), and r = —t. As a result we get the system 

Il9i f — ^'^xxx ^Zxxx ^‘^^Zxx T ^zZxx ^‘^'^xZx 48rc Wx 12z Wx T '2AwzZx T ^Zx , 

\ 37 = ^Wxxx — ^Zxxx — ^‘^WWxx + GzWxx — + QWxZx — 2AzWWx — 12w'^Zx + 2Az^Zx , 

for which we have the following proposition. 

Proposition 5. System (11911 is bi-Poisson 

~ + ‘^'^xZx — ‘^Zx — Qw'^Zx — Sz'^Wx — + Gw^z"^ — 2z‘^') dx 

= 1C5 f ^ {w^ — dx , 

with the compatible pair of Poisson structures 

D 0 


n = 


and K, having elements 


0 -D 


where 


ICn = (2Ai + 2D^ o As + o {D + iw) + (D - iw)-^ o (2Ai + 2A2D^ + 

T \2wxD Wx , 

^12 = {D — 4t(;) ^ o (l2ili + 2 B 2 D + 2B^D^ — AzD^ + 2D^^ + 12wxD~^Zx , 

/C 21 = -/Ct2 , 

IC 22 = (3Ci + 6L>2 o C 2 + iD^) o{D + 2z) ^ + {D - 2z) ^ o (3Ci + GC 2 D^ + 

+ l2zxD~^Zx , 

Ai = —Wxxx + Zxxx + 2AwWxx — ^2,WZxx “ ‘2^ZZxx + 16tC^ — Gzf + 24u;^Za; 

- IGGw'^Wx + 2‘iwzZx + ^z^Zx + 108u;^ - 2‘iw‘^z'^ — Az^ , 

A 2 = 4:Wx — 3za; — 12u;^ + 3z^ , 

Bl = —ZWxx — WxZx + 4:ZWWx , 

B 2 = SWxx — ^Zxx — 12wWx — ISzWx + ^ZZx + 36w‘^z — 4z^ , 
i?3 = 9wx — 6zx — 18w^ + 2z^ , 

(Si — Zxxx “b 8zZxx T (^^x 24^; Zx T 8z , 

<^2 = . 

It is obvious that PLis a Poisson structure. The operator /C is also a Poisson structure compatible 
with % but fractional decomposition of /C is a long expression to write explicitly. 

Transforming system ()17p by u = ^ and v = rs, we obtain 

w = rxxx + 3e(2r'^Sxx + ^rVxSx - srl - srrxx) + 3e‘^s‘^r‘^rx , 

— Sxxx T 8er(^sSxx T T 36 r s Sx , 

which is again a polynomial system. In this final form of system (|17|) . the coefficient of leading 
order terms is the identity matrix A = Ji and the parameter e is a factor of all nonlinear terms. 
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4. CONCLUSSION 

With the case of nondiagonal coefficient matrix of leading order linear terms, higher symme¬ 
try classification of (1, l)-homogeneous of weight 2 and 3 classes of two-component systems are 
completed. 

Although we naively completed our classification with nondiagonal matrix of leading order terms, 
all the systems obtained having a higher symmetry turned out to be systems admitting the Hopf- 
Cole transformation u = ln(t(;) 3 ; after which the coefficient matrix of leading order terms becomes 
the identity matrix Ji. By admission of a differential substitution, we mean that the system remains 
to be a system of differential (but not integro-differential) functions, after the substitution. 

On the other hand, the parameter e turned out to be essential for all systems obtained to be 
nontrivial: Although it is possible to set e = 0 in all systems obtained (e is nowhere in denominator), 
this reduction of all systems is triangular. 

Leaving aside the system m which is determined to be related to the well known Sasa-Satsuma 
system, we obtained three third-order symmetry integrable systems dUD-ClD and their master 
symmetries. These systems appear to be new. However, as C-integrable systems in the terminology 
of Calogero, these symmetry integrable systems are presumably useful nonlinear forms of linear 
systems. 

The second of the last three bi-Poisson systems (USD-dlZI) is related to a system given by Tsuchida. 
The remainin systems (|15l) and (jl7l) seem to be new. The Lenard-Magri schemes of these bi-Poisson 
systems are verified to work for the next two members of the symmetry hierarchy. Therefore these 
systems are integrable in the sense of infinite conservation laws. 
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